The probability distribution for non-Gaussianity estimators constructed from the 

CMB trispectrum 
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Considerable recent attention has focussed on the prospects to use the cosmic microwave back- 
ground (CMB) trispectrum to probe the physics of the early universe. Here we evaluate the proba- 
bility distribution function (PDF) for the standard estimator tvTi for the amplitude Tni of the CMB 
trispectrum both for the null-hypothesis (i.e., for Gaussian maps with Tni = 0) and for maps with a 
non- vanishing trispectrum {t^i 7^ 0). We find these PDFs to be highly non-Gaussian in both cases. 
We also evaluate the variance with which the trispectrum amplitude can be measured, <^Ar,ji), as 
a function of its underlying value, Tni. We find a strong dependence of this variance on Tni. We 
also find that the variance does not, given the highly non-Gaussian nature of the PDF, effectively 
characterize the distribution. Detailed knowledge of these PDFs will therefore be imperative in or- 
der to properly interpret the implications of any given trispectrum measurement. For example, if a 
CMB experiment with a maximum multipole of /max = 1500 (such as the Planck satellite) measures 
Tvii = then at the 95% confidence our calculations show that we can conclude Tni < 1005; as- 
suming a Gaussian PDF but with the correct Tni-dependent variance we would incorrectly conclude 
Tni < 4225; further neglecting the Tni-dependence in the variance we would incorrectly conclude 
Tni < 361. 

PACS numbers: 



I. INTRODUCTION 



Current observations of the cosmic microwave back- 
ground (CMB) and large-scale structure (LSS) provide 
a powerful probe of the physics of the early universe. 
As an example, the near scale-invariance of the primor- 
dial power-spectrum along with an upper limit to the 
inflationary gravitational-wave background can be used 
to rule out a few of the simplest models of inflation 
[1]. A measurement of the statistics of the primordial 
perturbations can provide an even more discriminating 
test of models of the early universe: all canonical single- 
field slow-roll inflation models predict that the perturba- 
tions are observationally indistinguishable from Gaussian 
[51[3]. Therefore any observed deviation from Gaussianity 
will rule out all canonical single-field slow-roll inflation 
models. However, non-canonical single- field models [1], 
multi-field models iSj, curvaton models [6], and models 
with sharp features [7 or wiggles [5] may produce larger 
departures from Gausianity. Measurement of the level of 
primordial non-Gaussianity has thus become one of the 
primary goals of CMB and LSS research. 

The majority of efforts to measure primordial non- 
Gaussianity from the CMB have relied on an estimator 
constructed from the bispectrum, the three-point corre- 
lation function in harmonic space [9j. However, most 
models that predict a significant level of non-Gaussianity 
also predict a non-zero connected trispectrum (the non- 
Gaussian part of the harmonic-space four-point function) 
[lOHTS] . and some efforts have been mounted to detect a 
primordial non-Gaussian signature from the trispectrum 



[TUl [T51 ITS] . In this way a constraint on the trispectrum 
amplitude provides unique information on a broad range 
of early-universe processes such as multi-field inflation 
models [T7|, the curvaton senario [TH], inflation models 
with non-standard kinetic terms |19j . and the influence 
of primordial cosmic strings |20j . 

The level of non-Gaussianity is often quantified using 
the 'local-model' through the non-Gaussianity parameter 
/ni defined by [S], 



$(.f) = 0(.f) + /„i [cf,{xf - (cj,')] 



(1) 



where ^(x) is the curvature potential and 4>{x) a Gaus- 
sian random field. Current limits from the CMB/LSS 
constrain the value to be |/ni| < 80 at 95% confidence 
level (c.l.) [2III23. The Planck satellite [3S] is expected 
to achieve a sensitivity of fn\ ^ 5. 

Constraints on the amplitude of the non-Gaussian 
local-model CMB bispectrum and trispectrum have very 
broad implications. Although various physical processes 
predict a range of values for f^i, it can be shown that all 
single-field models of inflation predict J4j 



/nl = ^("s- 1)> 



(2) 



where Ug is the slope of the primordial power spectrum. 
Current constraints to Ug [21] imply that all single-field 
models predict /„] ~ —0.008. Therefore, if the Planck 
satellite constrains /ni to be non-zero, we will be able to 
make the profound statement that all single-field models 
are disfavored by the data. A measurement of the ampli- 



tude of the local-model trispectrum^ , Tni, may lead to an 
additional test of the basic cosmological model. Recently 
Ref. [55] has shown that Tni respects the inequality 



Tnl > ^ifnl) , 



(3) 



independent of the underlying physics. A constraint on 
both /ni and t^i using the CMB may appear to violate 
Eq. ([3]) at the expense of actually violating translation 
invariance |26j . Therefore, a constraint on both f^i and 
t,ji provides a very broad test of some of the fundamental 
assumptions in our standard cosmological model. Given 
the wide-ranging impact of constraints on f^i and Tni, it 
is of great importance to report the significance of any 
constraint accurately. 

To date, most studies which use CMB observations to 
place constraints on /ni and Tni have used estimators that 
are constructed to have the minimum variance under the 
null hypothesis [TTJ [T^l [57] ■ In order to use these esti- 
mators to place meaningful constraints on /ni and Tni we 
must know the full shape of their probability distribu- 
tion functions (PDFs) as a function of /„!, Tni, and the 
maximum multipole l^ax of the observation. Thus, for 
example, we often evaluate or forecast the standard er- 
ror a with which a given measurement will recover the 
true value of /ni and Tni and then simply assume that 
the error is Gaussian. If so, then with crr„i = 100, for 
example, a measurement of TVa = 300 would represent a 
3a departure from Tni = 0, and a measurement tvj = 
would represent a 3cr departure fom Tni — 300. How- 
ever, if the PDF depends on the true value Tni, and if 
that distribution is non-Gaussian, then it may be that a 
measurement tvj — 300 could be easily consistent with a 
true value Tni = 0, while a measurement TVa = could 
be inconsistent with Tni = 300 with a confidence greater 
than "3tT." We will see below that something like this 
occurs with measurements of Tni. 

A calculation of these PDFs is particularly important 
for measurements of non-Gaussianity (as opposed, for ex- 
ample, for the CMB power spectrum), because the esti- 
mator is a sum over products of three (in the case of 
/ni) or four (in the case of t^) temperature measure- 
ments. This is unlike the power spectrum which sums 
over squares of temperature measurements. Suppose the 
temperature is measured in A^pix pixels. There are then 
~ TVpix terms in the /ni estimator and ~ N^^^ terms in 
the Tnl estimator (after restrictions imposed by statistical 
isotropy). While these terms may have zero covariance, 
they are not statistically independent; there is no way to 
construct N^^^ or N^^^ statistically independent quanti- 
ties from A^pix measurements. The conditions required 
for the validity of the central-limit theorem are therefore 



not met, and the estimators will not necessarily approach 
a Gaussian in the iVpix ^ 1 limit. 

Although previous studies have calculated the PDF in 
the case of /ni [2H1 [29] , the only property of the Tni es- 
timator that has been explored in the literature is how 
the variance for the null-case scales with the maximum 
observed multipole, ^niax [I2[I11[IS]. In order to address 
how well CMB observations can estimate Tni we calculate 
the PDF i-'[T\n;Tni, ^jnax] — the probability that a given 
measurement with resolution Imax will return a value tv^ 
given that the underlying theory has a value Tni — using 
numerous Monte Carlo realizations of an ideal (no-noise) 
flat-sky map in the Sachs- Wolfe approximation. In order 
to both generate the maps and apply the estimator to 
them we use a fast-Fourier transform (FFT) algorithm 
described in Appendix A of Ref. [29]. Lessons learned 
about P[TnhT^i,li-na.x] ill this ideal case help to interpret 
and understand current /forthcoming results and assess 
the validity of full-experiment simulations. 

Our simulations show that P[TVd; Tni, Zmax] is highly 
non-Gaussian for all values of Tni, including the null case. 
Additionally, our simulations allow us to derive, for the 
flrst time, how the variance of this distribution depends 
on the underlying value of Tni. Neglecting this depen- 
dence, Ref. [14] concludes that the signal-to-noise^ of 
this estimator appears to scale as Tni^, becoming more 
sensitive to non-Gaussianity for large Tni than an esti- 
mator using the CMB bispectrum. We show that the 
dependence of the variance on the underlying value of 
Tnl is significant, finding that the sensitivity of this es- 
timator to local-model non-Gaussianity is always weaker 
than that of the bispectrum estimator, and it approaches 
a constant for Tni ?J lO^/^max- 

Knowledge of both the variance and shape of 
P[fni',Tni,lma.x] is ncccssary to assign proper confidence 
levels (c.l.) to constraints. For example, if the Planck 
satellite (^max — 1500) measures tvj = then at the 
95% c.l. our calculations show that we can conclude 
Tnl < 1005. If we assumed a Gaussian PDF but with 
the correct Tni-dependent variance we would incorrectly 
conclude Tni < 4225. If we also neglected to include 
the Tni-dependent variance we would incorrectly conclude 
Tnl < 361. 

This paper is organized as follows. In Sec. [llj we dis- 
cuss how to construct the minimum-variance estimator 
T^ using the CMB trispectrum under the null-hypothesis. 



In Sec. Ill we apply this estimator to the local-model for 



^ The local-model trispectrum can be defined by using Eq. M with 



the identification Tni = /nl 



non-Gaussianity. In Sec. [Ill A| we present our results for 
the PDF in the null (Tni = 0) case. Sec. Ill B presents our 
results for the non-null case and gives a fitting formula 



^ The signal-to- noise is defined to be S/N = T^^i/ar^_^^. In the case 
of Gaussian noise with a variance independent of r^i this is a 
measure of the fractional error in a constraint to r^i. However, 
in the case of t^, the noise is neither Gaussian nor independent 
of Tnl so that the quantity Tni/o-^^i is only an approximation to 
the significance of a constraint on r^i. 



for P[t,ii; t^i, ^max]- In Sec. IV we summarize our results 



B. The minimum-variance trispectrum estimator 



II. NON-GAUSSIANITY ESTIMATORS 

CONSTRUCTED FROM THE CMB 

TRISPECTRUM 

A. Formalism 

We assume a fiat sky to avoid the complications 
(e.g., spherical harmonics, Clebsch-Gordan coefficients, 
Wigner 3j and 6j symbols, etc.) associated with a spheri- 
cal sky, and we further assume the Sachs- Wolfe limit. We 
denote the fractional temperature perturbation at posi- 
tion 6* on a flat sky by T{6) and refer to it hereafter simply 
as the temperature. 

The field T{d) has a power spectrum Ci given by 



TrTj- ) — ilSf . r „Ci, 

/i h / li+l2,0 " 



(4) 



Each distinct quadrilateral ^i + /2 + ^3 + ^4 = gives 
an estimator for the trispectrum with some variance. 
Adding the individual estimators with inverse-variance 
weighting gives the minimum-variance estimator. 



^1 = f^T,0 X! 

li+l2+h+U=0 

X T{1'iX,i'3X)-<Jt.q{T) 



Tf Tf Tf Tf 

1-1 t2 t3 U 



m^Ci,Ci,Ci,Ci, 



(10) 



where we have subtracted off the unconnected (Gaussian) 
part of the trispectrum, {T)q, and the inverse variance 
is given by 



'T,0 



E 

il+i2-|-i3-|-'4 



T{li,l2,hih) 



m-^Ci,Ci,Ci,Ci, 



(11) 



where 2 < j^^j < ^j, 



where 17 = 47r/sky is the survey area (in steradians). 



Tr= I d^ee-''-'^T{e) 



A^n 



J2e-'''m, (5) 



is the Fourier transform of T{9), and ^f , j: g is a Kro- 

necker delta that sets h = — ^2- In the limit of small de- 
partures from Gaussianity, C/ is also the power spectrum 
for T{9), which for a scale-invariant primordial power 
spectrum with amplitude A (A ~ 10~^") is given by 



Ci = 



2t:A 



The trispectrum is defined by [TTl [T^ 



(6) 



and for the local-model, 



nWMMM) = P!±{\h + h\) 



+ p}:!:(.\h + h\) + p!±i\h + Q), 



(8) 



where 



P!±{\h + r2|) ^ 4C|,- ^,- I [Ci,Q, + Q,Q, 



(9) 



III. THE PDF OF ^i FOR THE LOCAL-MODEL 

We now restrict our attention to the local family of 
non-Gaussian models [see Eq. (II])] in which the temper- 
ature T{6) has a non-Gaussian component, 

ne) = m + ?../^i[mf {mY)], (12) 

where we have chosen the normalization Tni to correspond 
to amplitude of the non-Gaussian part of the Newtonian 
potential four-point function'^. To zeroth order in Tni, 
the power spectrum and correlation function for T{9) are 
the same as those for t{6). Note that T{9) is, strictly 
speaking, the temperature fluctuation, so (T{9)j = = 
T- 

The temperature Fourier coefficients can be written 
Tf^tf+V^^iS with 



5tj ^ ^ 



E 

ll+l2=l 



ti h 



(13) 



Formally, the sum goes from 1 < |^i| < oo, but for a 
finite-resolution map, the sum is truncated at some /max 
such that the number of Fourier modes equals the number 
of data points. 



Applying the estimator in Eq. ( 10 1 to the local-model 
trispectrum [Eq. Q], we obtain 



?^ = 2 Y. 

l<|L|<2i„ 



Cl 



E 



'l '2 



h+l2+L = 



n^Ci, 



4,0 (T) 6^4) 



Due to statistical isotropy, the trispectrum is nonvan- 
ishing only for ^i -f- ^2 + ^3 -f ^4 = 0, that is, only for 
quadrilaterals in Fourier space. 



^ This differs by a factor of 5/6 the usual definition which is given 
in terms of the Bardeen potential 1301 . 



where we can now write 



{T)a = 24,0 E ^^ E 



|L|>1 li+h+L=0 



(15) 



A. The PDF of Tni under the null-hypothesis, Tni = 

Under the null hypothesis (thi = 0) we apply t^ to a 
purely Gaussian CMB temperature map. As shown in 




max 



FIG. 1: The scaling of the variance of t^ under the null hy- 
pothesis Tni = 0. Each red point shows the result of a Monte 
Carlo simulation for 1000 realizations. The black line is a fit 
to those simulations given by Eq. \1G\ . 

Fig. [1] our Monte Carlo simulations find that the vari- 
ance (To of this estimator under the null hypothesis as a 
function of the maximum multipole ^max included in the 
analysis is well-fit by a power-law 



<^OV^riiaxj — 



1.74 X 10" 

A21 4 



(16) 



This scaling compares well with the results of previous 
work [ni[Tl]. 

The simulations also allow us to calculate the full shape 
of the PDF of this estimator under the null hypothesis. 
Since the number iVpjx of measurements of the CMB tem- 
perature map is much less than the number (~ -^pix) of 
terms in this estimator, the standard central-limit theo- 
rem does not apply so that the PDF P[fvj; Tni = 0, /max], 
is not necessarily Gaussian in the iVpix ^ 1 limit. The 
simulations demonstrate that the PDF is, in fact, highly 
non-Gaussian as shown in Fig. [5] The asymmetry of the 
PDF can be explained by referring to the expression for 
the estimator in Eq. (14). There it can be seen that the 



above. By calculating the PDF, P[rni;Tni — 0, /max], for 
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FIG. 2: The black-solid curve shows the PDF of fJTi under the 
null-hypothesis, Tni = for ?niax = 50 and calculated with 10^ 
realizations. The red-dashed curve shows a Gaussian PDF 
with the same variance. The upper panel shows the PDF on 
a linear scale, the low panel on a logarithmic scale. When 
scaled by its variance, the PDF is identical for /max = 100 
showing that it takes on a universal form in the /max 3> 1 
limit. We give a fitting formula for this PDF in Eq. ( 17l. 



several values of /max we find that when scaled by its vari- 
ance it takes on the universal shape shown by the black 
curve in Fig. [2] The PDF is well fit by the formula, 



P^V^TrA = 0,/max] = 



1 



(17) 



, 1 ^nl/" 



3 "^ ^ T„i/cr>a;p, 



where a is the variance of the estimator, given in Eq. ( 16 1 , 
the normalization N is given by 



N = 2apT 



n + l 






(18) 



estimator is bounded from below but unbounded from 



where P is the Euler Gamma function, and Ki is the 
modified Bessel function of the first kind. We find that 



the PDF is best fit by the parameters n — 3, Xp — —0.13, 
CTo = 0.64, and c = 0.' 



B. The PDF of f^i with Tni / 

For Tni ^ the non-Gaussianity in the CMB 
map imparts further non-Gaussianity to the shape of 
^[tVli; ■''ni, ^max]- In addition to this, the variance has a 
strong dependence on Tni so that when t^i and ^max are 
large enough the ratio T-^\/(JTan which approximates the 
S/N of the estimator, approaches a constant value. 

In order to investigate how the variance of TVa depends 
on Tni it is useful to expand it in powers of Tni. Given 
that T; is linear in Tni and that t^\ is quartic in T; the 
expansion includes terms up to Tni^: 

^ = To + V^Tl + TnlTa + Tj'^Ti + Tnl2r4, (19) 

where each % ^ St^^ . We give explicit expressions for 
the 7i in Appendix A. 

Since only cross-correlations which include even prod- 
ucts of t are non-zero, the variance of tvj is given by 



^[ATnl]^^ = Cro+Tnl(CT?+ 0-^,2) 



(20) 



Tni {(^2 + Cro,4 + <^1,3) 

^ 3/^2 .2 \ ,^ A 2 
T^nl (.CTg + (72, 4J + ^nl 0^4 



where, for example, CTq 2 denotes the covariance between 
% and T2 and At^ = tvS - (tI^). 



Calculating the terms that appear in Eq. (21) using 
Monte Carlo simulations we find that for Tni < lO'* and 
^max < 10'* only a\ and a\ significantly contribute. The 
variance is then well approximated by 



[A?^ 



2 

T'nlfi 



2 2 
Tni 0-2; 



where CTq is given in Eq. ( 16 1 and 



0.028 
4/ 2 
0.23. 



(21) 

(22) 
(23) 



The results of our Monte Carlo simulations are shown in 
Fig.i 



The scaling given in Eq. (21) shows that for a large 
enough Zniax the variance of the estimator scales as Tni^ 
so that the ratio TVLi/ff^j becomes constant for Tni ^ 
0.1/(AZniax^)- A similar scaling is observed with the 
minimum-variance null-hypothesis estimator using the 
CMB bispectrum [28l [29] . Neglecting the dependence of 
the variance on Tni, previous work [Tl] claimed that for 
large enough Tni the minimum-variance null-hypothesis 
estimator using the CMB trispectrum would be more 
sensitive to a local-model non-Gaussian signal than an 
estimator using the CMB bispectrum with the relation- 
ship Tni = /ni • Given the dependence of the variance on 




max 



FIG. 3: The variances a\ and a\ as a fraction of the zeroth- 
order variance g%. The points are the results of the Monte 
Carlo si mul ations and the curves show the power- law fits given 
in Eqs. pSl) and (l23l). 




20 40 60 80 100 



/nl (= \/^) 



FIG. 4: The ratio fni^/(^j^i using the CMB bispectrum (black, 
from Ref. [29j ) and Tni^/a^^^ using the CMB trispectrum 
[blue, Eq. ([2l|] as a function of Tni = /ni^. These ratios can 
be interpreted as an estimate for the S/N for a constraint to 
local- model non-Gaussianity in the CMB. The dashed curves 
show the scaling of the S/N without taking into account 
the dependence of the variance on /ni and Tni; the solid 
curves show the correct S/N scaling. As in Ref. [14], from 
the dashed curves we would (incorrectly) conclude that the 
trispectrum estimator is more sensitive to a non-Gaussian sig- 
nal for ^TVu > 40. 



Tni our calculations demonstrate that this is not the case, 
as shown by the solid curves in Fig. |4J 

When tJ^ is applied to a map with t^i ^ then 
the non-Gaussianity in the map imparts additional non- 
Gaussianity to P[T^i;Tni,lnia.x]- We are interested in cal- 
culating the shape of the PDF for an experiment such 
as Planck which has /max — 1500. Although, in prin- 
ciple, it is possible calculate the PDF for large /max, it 
is computationally demanding especially given the large 
number of realizations we must generate in order to ex- 
plore the tails of the distribution. The computation can 
be simplified since we find the PDF is 'self-similar' in 
the sense that its shape depends on the ratios (Jg/al and 



crg/cr|. Given Eqs. (|16|), (|22|), and (|23| this implies that 



max' 



the shape depends on the combination t^iI^ 

Using this fact it is straightforward to calculate the 
PDF for a moderate value of /max (we used /max = 50) 
and then scale the PDF to a larger value for various 
choices of t^i. We find that P[T\a': TniJmax] is well fit 
by the formula used to fit P[t^;t„i — 0, /max], given by 
Eq. ( |17[ ), with parameters n, ap, Xp, and c which now de- 
pend on /„] and a is the variance of the estimator given by 
Eq. ( [21] ). In Fig.lSlwe show how these parameters depend 
on Tni for /max = 600 (dotted), /max = 1500 (dashed) and 
^max ~ 3000 (solid). We find that as Tni increases the 
asymmetry of the PDF increases with the power-law in- 
dex of the PDF for t\j < Tni growing from n = 3 to n = 4 
for /max = 1500 and n = 4.24 for /max = 3000. 

Our knowledge of the full shape of P[fni; Tni, ^max] now 
allows us to properly assign confidence levels (c.L). If a 
given CMB observation with /max yields a value of (T\j)obs 
we assign the 95% c.l. by finding the value of (Tni)±c.i. 
which satisfies the integral equation 



0.95 = 



(^i)± 



(f;ri)„b 



^[■^ni; (Tnl)±c.l., 'maxl^Tnl 



(24) 



where (jni)± is the solution to the equation 



-P[('''nl)±; ('''nl)±c.l., ^max] — ^[('''nOobs; ('''nl)±c.l. , ^max] 



If, for example, an experiment with /max = 3000 mea- 
sures (^) = 100 then at the 95% c.L, Tni = lOOtgf- 
If we assumed a Gaussian PDF with the Tni-dependent 



Tni 



100+3180. 



variance given by Eq. (21 1 we would incorrectly conclude 
Finall 



y, if we also neglected to include 
the Tni-dependent variance we would incorrectly conclude 
Tni = 100 ± 92. 



IV. DISCUSSION 

In this paper we have shown that the PDF for non- 
Gaussianity estimators using the CMB trispectrum can- 
not be assumed to be Gaussian, since the number, r^ 
TVpj^, of terms used to construct these estimators greatly 
exceeds the number A^pix of measurements. The 99.6% 
confidence-level interval cannot safely be assumed to be 




10 20 30 40 50 



v^ 



1 



FIG. 5: The scaling of the parameters for the fitting formula 
in Eq. (17 1 with y/Tni for /max = 600 (dotted curves) /max = 
1500 (dashed curves) and /max = 3000 (solid curves) . For 
Tni 7^ the non-Gaussianity of the map imparts additional 
non-Gaussianity to the PDF. As Tni increases the asymmetry 
of the PDF increases with the power-law index of the PDF for 
Tni < Tni going froui n = 3 to n = 4 in the case of /max = 1500, 
and n = 4.24 in the case of /max = 3000. 



3 times the 68.2% confidence-level interval. We found 
that the PDF of standard minimum-variance estimator 
Tv^ using the CMB trispectrum constructed under the 
null hypothesis is well- approximated by a distribution 
given by Eqs. (171 and (21). This distribution is expo- 



nentially suppressed for va 



ues Tni < Tni and enhanced 



for values f^\ > Tni, relative to a Gaussian with the same 
variance. We calculated how the parameters of this fit- 
ting formula depend on Tni for /max = 600, /max = 1500 
and /max = 3000 as shown in Fig.[5j We also find that the 



(25) non-Gaussianity of -P[tVli; Tni, ^max] is greater for Tni 7^ 0. 



We have calculated, for the first time, how the variance 
of T\j depends on the underlying value of Tni, as shown in 
Eq. (21 1 . Previous work neglected this dependence lead- 



ing to the incorrect conclusion that for large enough Tni 
and /max a non-Gaussianity estimator constructed from 
the CMB trispectrum would have a larger S/N than an 
estimator constructed from the CMB bispectrum [14] . 
When the Tni dependence in the variance is included the 
S/N of the estimator constructed from the CMB trispec- 
trum becomes constant for Tni > 0-l/(^'max)- ^s a re- 
sult, the estimator constructed from the CMB bispec- 
trum always produces a larger S/N, as shown by the 
solid curves in Fig. [4J 

These results have important consequences for future 
constraints to Tni measured from the CMB. As discussed 
in the Introduction, future constraints to both /ni and Tni 
may imply wide-ranging conclusions about the physics of 



the early universe. In particular, a non-zero measurement 
of /ni would rule out all single-field inflation models and a 
constraint to t^i probes basic physical assumptions of the 
early universe, such as translation invariance, by testing 
the consistency relation, Tni > /ni /2 around the surface 
of last scattering [Ml |26] . If we suppose an experiment 

with Zmax — 3000 measures /„] — 20 and t^ = then at 
the 95% c.l. our calculations show that we can conclude 
Till < 200, violating the consistency relation Tni > /ni /2 
at the 95% c.l. If we assumed a Gaussian PDF for fn\ but 
with the correct Tni-dependent variance found in Eq. ( 21 ) 



we would incorrectly find Tni < 1000, leading to the false 
conclusion that Tni > /ni /2 is consistent with the data. 
If we also neglected to include the T,ii-dependent variance 
we would incorrectly find Tni < 90. 

The non-Gaussian shape of the PDF of fn\ is im- 
portant even for current constraints. The current pub- 
lished constraints on Tni from the Wilkinson Microwave 
Anisotropy Probe (WMAP) [H [21] (^max = 600) are 
Tni/lO^ = 1.68 ± 1.31 at the 68% c.l. [IS]. The error 
quoted in this constraint is estimated without taking 
into account the full shape of the PDF. Although our 
results are not directly applicable to this case since they 
do not take into account several details of the WMAP 
analysis (such as a full CMB transfer function and the 
noise properties of the observations) they do allow us 
to discuss how using the correct PDF would qualita- 
tively change the confidence levels. For Z^ax = 600 our 
calculations show that when assuming a Gaussian PDF 
along with a Tni independent variance the constraint is 
Tni/lO"* = 1 ± 0.1; the full PDF shows the actual con- 
straint to be Tni/10'' = llo;2^4)^g^. This indicates that a 
complete treatment of the confidence levels for the cur- 
rent constraint on Tni is both asymmetric and has a larger 
range than the constraint that is quoted in Ref. [IB]. 

The results presented here are made within the flat- 
sky, Sachs- Wolfe approximation. As such our conclusions 
should be taken as an order-of-magnitude estimate of 
P[fni', Tn\, ^niax] Calculated on the full sky and with the full 
transfer function (see Ref. [28] for a further discussion). 
However, we note that a comparison between the exact 
and approximate scaling of the S/N with Zniax shows the 
agreement to be better than an order of magnitude |27j . 

In this paper we have concentrated solely on the non- 
Gaussian local-model, defined in Eq. ([l]). Although the 
quantitative results will differ for other models of non- 
Gaussianity, we expect that the qualitative conclusions 
will remain unchanged. The non-Gaussian PDF for t^i 
results from a breakdown of the central-limit theorem 
due to the large number of terms used in the estimator 
compared to the number of independent measurements. 
Therefore, the work presented here shows that estimators 
for the CMB trispectrum amplitude cannot be assumed 
to have a Gaussian PDF. Rather, one must carefully ex- 
plore the full shape of the PDF before assigning the sig- 



nificance of any particular measurement. Similar consid- 
erations may also need to be considered, for example, in 
measurements of things like weak gravitaitonal lensing, 
departures from statistical isotropy [32^, and the like, as 
the magntiudes of many of these effects are determined 
in practice by the trispectrum. 
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Appendix A: The expansion of Tni in Tni 

Here we write down the explicit formulas for the ex- 
pansion of T\a, the minimum- variance estimator using the 
CMB trispectrum constructed under the null hypothesis, 
as written out schematically in Eq. (19). The standard 



estimator applied to the non-Gaussian local-model, de- 
fined by Eq. pj, can be written 



(Al) 



Not: 



Tni 



r^i + (t'^.o 01 G = 










Tj-Tj- 

=0 


ing that Tf = if + ^^T^iStj. we have 





2imax 



|L| = -2/„„x, L#0 
Tnl(|A2p+AiA;+^3At) 



+ tJ'^ {A2AI + A:iAl) + Tni^lAart - (?^) 
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